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We investigate the gravitational lens effect caused by a non-spherically distorted compact 
object. The non-spherical property of the gravitational potential is modeled by a quadrupole 
moment. Under the assumption that the quadrupole contribution is small, we solve pertur- 
batively the lens equation and obtain the image positions and the amplification factors. We 
show that the separation angle of two major images is only slightly changed by the existence 
of the quadrupole contribution, whereas the difference of the amplification factors may be 
significantly modified. Our results indicate that even a tiny non-spherical distortion of the 
lens potential may cause significant amount of flux anomalies in the lensed images. 

§1. Introduction 

The gravitational lens is an important tool in astrophysics for probing mass dis- 
tributions and determining the cosmological parameters.^ Although previous stud- 
ies of gravitational lensing have been mostly based on simple, spherically symmetric 
lens models, some generalizations which include rotation® - ® and higher-order gen- 
eral relativistic effects^ >EJ have also been developed. In this paper, we investigate 
another generalization, i.e., the effect of non-spherical distortion of the gravitational 
potential in a compact lens object on the gravitational lensing. 

A pioneering work on the non-spherically deformed compact lens was done 
by AsadaP He studied analytically the gravitational lensing caused by a non- 
spherically deformed star, where the non-spherical property of the gravitational po- 
tential was modeled by a quadrupole moment. He employed a rigorous analytic 
approach and presented solutions of the image positions for a source on the principal 
axis, an expression of the caustics and the critical curves. However, because of the 
highly non-linear nature of the lens equation, general solutions of the image positions 
for a off-axis source and their amplification factors have not been presented. 

In this paper, we take a more practical approach, namely the perturbative ap- 
proximation, and solve analytically the lens equation and obtain the approximate 
solutions up to the lowest order of the quadrupole moment. We also calculate the 
Jacobian of the lens mapping and indicate how the quadrupole moment changes the 
amplification factors. 
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§2. Lens equation 

In this section, we summarize the lens equation for a compact object with a 
quadrupole moment. See also Asada® for detail. The gravitational potential of a 
non-spherically deformed compact object is modeled by use of a monopole and a 
quadrupole moment: 

i , , GM G ( iV T , . 

4> = ^ + 4>2 = — - (^3— - -g-J iij, (2-1) 

where 



M = y pcTx, 7y = y pxiXjd^x, (2-2) 

and /O is the mass density of the lens object. Using the above form of the potential, 
the deflection angle 2-vector a. is calculated as follows: 



AGM f 8G / ^'^f g) , 



where the 2-dimensional vector £ denotes the image position, and 

Qij = J p hCiXj - ^SiAxA d 3 X (24) 

denotes the trace- free quadrupole moment. Without loss of generality, we can assume 
the normalized quadrupole moment is diagonalized as follows, 

® lj ^ 2GM*D L D LS Ql3 = ( -e ) ' (2 ' 5) 

where D$, Dl, Dls are the angular diameter distances from the observer to the 
source, from the observer to the lens, and from the lens to the source, respectively. 
Hereafter, we assume e > 0. Finally, the lens equation is 

- x X - e (x2 ~ 3y2)x (2-6) 

Px ~ x 2 + y 2 (x 2 + y 2 )3 ' 1 J 

r - y - ( 3x2 - v 2 )v (0 7 x 

Pv-y x 2 + y 2 e ( x 2 +y2) 3> V) 

where (3 = (f3 x ,(3 y ) and 6 = (x,y) are the source and the image positions, respec- 
tively, normalized by the Einstein radius 9e = V 4GM D is I 'c 2 'D lD s ■ 

§3. Perturbative approach 

A comprehensive, algebraic study of the above lens equation Eqs. (|2-6p and 
(|2-7p was done by AsadsP using polar coordinates. He showed the lens equation for 
this type as a single real lOth-order algebraic equation and gave analytic solutions 
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for a source located exactly on the principal axes. He also estimated the order of 
magnitude of the normalized quadrupole moment as 

^WffsW « W^V( '- ,)', (3-1) 
VMyVl0 6 kmy V / VlOkms" 1 / V ' 

where R denotes the typical size of a lens object, Re is the Einstein ring radius, and 
v is a rational surface velocity on the equatorial plane of the lens star. Then, for a 
nearby solar- type star at 10 pc, i.e., R ~ 10 6 km, Re ~ 10 7 km, with a rotational 
velocity of v ~ 10 km s _1 which is much faster than that of the Sun, the estimated 
value is e ~ 10~ 5 , which is sufficiently smaller than unity. 

Because of the higher polynomial nature of the lens equations, however, general 
analytic solutions for a source not on the axes have not been given yet (or cannot be 
given for more than 4 images). Also, the physical properties of the lensed images, 
such as the image amplification factors, are not yet investigated. 

Instead of such an exact algebraic treatment, in this paper, we take another, a 
more tractable approach to the quadrupole lens equations. Using the fact 0<e<l, 
we employ a perturbative approach to the lens equations and obtain approximate 
solutions to the image positions (x, y), up to the lowest order of the eigenvalue e of 
the quadrupole moment. 

Concerning the perturbative approach to the gravitational lensing, for exam- 
ple, AlarcP™ wrote a series of papers in a slightly different context from ours. 
His main interest is the non-spherical perturbations of extended lens models and 
elongated arc images, whereas in this paper we concentrate our attention on a non- 
spherically distorted compact lens model and multiple images. 

We start from the following set of the lens equation, which is obtained by clearing 
the fraction of Eqs. (jM)) and (J2ZTJ): 

(x 2 + y 2 ) 2 {(x 2 + y 2 )(p x - x ) + x} = -e(x 2 - 3y 2 )x, (3-2) 

(x 2 + y 2 ) 2 {(x 2 + y 2 )(p y -y)+y} = -e(3x 2 - y 2 )y. (3-3) 
Setting e = 0, we obtain the zeroth-order solutions, 

x = x± = f ± p x , y = y^f ± (3 y , (3-4) 

where 



(3-5) 

and (3 = *//3j~+~/3| . The above solutions are just the solutions for a point mass lens. 

We may also find a trivial solution (x,y) = (0,0), but it is inadequate because the 
denominators in the original lens equation Eqs. (|2-6p and (|2-7p vanish. 
Next, we put the following form into Eqs. (|3-2p and (|3-3j) . 

x = x +xi, 2/ = 2/o + 2/1, (3-6) 
then, up to the linear order of x\ and y\, we obtain 

{xl - 3>ypx Q 
(4 + y 2 f 



{2x QXl + 2y oyi )(p x - x ) - {xl + yl) Xl +x x = -e \°_ 2 , y °' , (3-7) 
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(2x0*1 + 2y oyi )(f3 y - yo) - (x 2 + yl) Vl + yi = -e ^*° ■ (3-8) 

The solutions are expressed in terms of /3 X and j3 y as follows: 



*± - *± + x ± - f ±0 + c ( 4 ^- 3 / 32 )(/ ± ) 2 - 1 » (3 . 9 ) 
+ + 4. + (30 2 -4#?)(/±) 2 + l 

= vt + = + vW + ix^ + a) ^- (340) 

In the case of a point mass lens model (e = 0), the number of the images is two. 
As shown by Asadaj*^ the number of the images by a quadrupole lens (e / 0), which 
depends on the value of e and the position of the source, is generally more than two, 
in most cases four. We can obtain another, new solutions for e > case, which are 
not existent for a point mass lens CcLS6 ; clS cl perturbation around the trivial solution 
(x,y) = (0,0), namely, 

x = + xx, y = + yi. (3-11) 
Up to the lowest non-trivial order, the lens equation is 

(x 2 + y 2 ) 2 xi = -e(x 2 -3y 2 )xi, (3-12) 

{xl + ylfy^-e^xl-yDyL (3-13) 

In the case of x\ ^ 0,yi ^ 0, Eqs. (|3-12|) and (|3-13|) imply x\ = — y 2 , which does 
not have a real root. In the case x\ ^ 0, y\ = 0, Eq. (|3T2p implies x\ = — e, which 
again does not have a real root since we have assumed e > 0. Finally in the case 
x\ = 0, y\ ^ 0, Eq. (|3-13|) reduces to y\ = e, from which we obtain 

x = xi = 0, y = yi = ±Ve- (3-14) 

The above solutions are always on the y-axis, independent of the source position 
(PxiPy), and order of 0{y/e), not 0(e). Therefore, we perform one more iteration 
and obtain slightly higher order solutions in the following form, 

x = + x 2 , y = + yi + y 2 = ±yfe + y 2 . (3-15) 

Up to the lowest order of x 2 and y 2 , the lens equation is now 

ViPx + yfx2 = 3ey 2 x 2 , (3-16) 

yt^y-yi) + 5yty 2 = 3ey 2 iy2 , (3-17) 

and the solution is 

x = x 2 = ~eP x , y = yi+y 2 = ±y/e(l + ^e\ -~eP y . (3-18) 

These images always appear very close to the y-axis, are very dim because the 
amplitude of the amplification factor is 0(e 2 ), and disappear when e — > 0. We call 
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Table I. Comparison of our approximate solutions with numerical ones for B x = 



(- 


Bx 


Hy 


"■'Hum 




1 1 ) 7/*,.,™ 
V / i/num 


(2) £/appr 


1(2) - fill 


0.01 


0.0 


0.2 


0.0 


0.0 


1.10087 


1.10091 


3.9 x 10~ b 








0.0 


0.0 


-0.89881 


-0.89891 


1.0 x 19~ 4 








0.0 


0.0 


09950 


099500 


3.7 x 10 -6 








0.0 


0.0 


-0.10157 


-0.10150 


6.8 x 10~ 5 


0.01 


0.0 


0.5 


0.0 


0.0 


1.27780 


1.27782 


1.8 x 10 _b 








0.0 


0.0 


-0.77262 


-0.77282 


2.0 x 10" 4 








0.0 


0.0 


0.09809 


0.09800 


8.5 x 10~ 5 








0.0 


0.0 


-0.10327 


-0.10300 


2.7 x 10" 4 


0.02 


0.0 


0.2 


0.0 


0.0 


1.09668 


1.09684 


1.6 x 10" 4 








0.0 


0.0 


-0.89241 


-0.89284 


4.2 x 10 -4 








0.0 


0.0 


0.14084 


0.14084 


2.5 x 10" 7 








0.0 


0.0 


-0.14510 


-0.14484 


2.7 x 10" 4 


0.02 


0.0 


0.5 


0.0 


0.0 


1.27479 


1.27486 


7.3 x 10" 5 








0.0 


0.0 


-0.76402 


-0.76486 


8.4 x 10" 4 








0.0 


0.0 


0.13802 


0.13784 


1.8 x 10" 4 








0.0 


0.0 


-0.14878 


-0.14784 


9.5 x 10" 4 



them as the "minor" images, and may safely neglect them in the analysis of the 
image amplifications. 

In order to check the accuracy of our approximate solutions Eqs. (|3-9p . (|3-10p . 
and (|3-18p . we compare them with numerical solutions. Table HI and HI1 show the result 
for some sets of typical values of the parameters. Table [I] shows the calculated values 
for the source on the y-axis, i.e., /3 X = 0, and Table HI1 for the source on the x-axis, 
/3 y = 0. The maximal error is about 10 -4 , which is the same order of magnitude as 
0(e 2 ). 



Table II. Comparison of our approximate solutions with numerical ones for B y = 



e 


B x 


Pv 


(l) ^num 


(2) X&ppr 


1(2) -(1 


)l 


(3) |/num 


(4) J/appr 


1(4) - (3)| 


0.01 


0.2 


0.0 


1.10902 


1.10906 


3.8 x 10" 


-■ 


0.0 


0.0 


0.0 








-0.91097 


-0.91106 


9.7 x 10" 


■5 


0.0 


0.0 


0.0 








0.00102 


0.00100 


2.0 x 10" 


■5 


0.10048 


0.10050 


1.7 x 10 -5 








0.00102 


0.00100 


2.0 x 10" 


■5 


-0.10048 


-0.10050 


1.7 x 10" 5 


0.01 


0.5 


0.0 


1.28372 


1.28373 


1.8 x 10" 


■5 


0.0 


0.0 


0.0 








-0.78855 


-0.78873 


1.9 x 10" 


■4 


0.0 


0.0 


0.0 








0.00254 


0.00250 


3.8 x 10" 


■5 


0.10035 


0.10050 


1.5 x 10" 4 








0.00254 


0.00250 


3.8 x 10" 


■5 


-0.10035 


-0.10050 


1.5 x 10" 4 


0.02 


0.2 


0.0 


1.11299 


1.11314 


1.5 x 10" 


-4 


0.0 


0.0 


0.0 








-0.91676 


-0.91714 


3.8 x 10" 


■4 


0.0 


0.0 


0.0 








0.00208 


0.00200 


8.1 x 10" 


■5 


0.14281 


0.14284 


2.5 x 10 -5 








0.00208 


0.00200 


8.1 x 10" 


■5 


-0.14281 


-0.14284 


2.5 x 10" 5 


0.02 


0.5 


0.0 


1.28662 


1.28669 


7.0 x 10" 


■5 


0.0 


0.0 


0.0 








-0.79598 


-0.79669 


7.1 x 10" 


■4 


0.0 


0.0 


0.0 








0.00516 


0.00500 


1.6 x 10" 


■4 


0.14241 


0.14284 


4.2 x 10" 4 








0.00516 


0.00500 


1.6 x 10" 


■4 


-0.14241 


-0.14284 


4.2 x 10 -4 
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§4. Amplification factor 



Once we express the solutions for the image positions (x, y) in terms of the 
source position (f3 x ,f3 y ), we can directly calculate the amplification factor from the 
Jacobian 

d(x ± ,y ± ) 



A 



± 



det- 



(4-1) 



d{p x ,Py) 

where A + and A~ represent the amplification factors for the images with positive and 
negative parity, respectively. The calculation is straightforward, and from Eqs. (|3-9p 
and (|3-10p we obtain the following results up to the linear order of e: 



A + = Ai 



m ~ Hi) W 2 + 6) - W 2 + 4) V / ^+4 



(4-2) 



where 



A' 



A- n 



4(/3| - /3 2 ) /3(/3 2 + 6) + (/3 2 + 4) v^TI 



/33(/3 2 + 4) 



A ± 



f3 2 + 4±(3 



4/3^^2 + 4 



(4-3) 



(4-4) 



represent the amplification factors for a point mass lens, namely e = case. 

We can also define the Pade approximants for A ± of order [0/10 with respect 
to e as follows: 



At 



1 + e- 



m - p 2 y ) m 2 + 6) - (/3 2 + 4) 



/? 3 (/3 2 + 4) K//32 + 4 + /3 



(4-5) 



a: 



A- n 



1 + e- 



4(/3 2 - /3 2 ) /3(/? 2 + 6) + (/3 2 + 4)v^T4, 



V 



/3 3 (/3 2 + 4) K//3 2 + 4-/3 



(4-6) 



As shown in Table UTIl the Pade approximants Eqs. (|4-5p and (|4-6p generally give 
the approximate amplifications with better accuracy of about 1% or less. Hereafter, 
we use Eqs. (|4-5p and (|4-6p as our approximate formulae for the amplification factors. 



*' Generally speaking, the Pade approximant of order [m/n] is an approximation of a function 
by a ratio of two polynomials of order m and n. Then, for a function A(x) whose Taylor expansion 
is given by A(x) ~ Ao(l + a\x +■■■), the Pade approximant of order [0/1] is simply Aq/(1 — a\x). 
At the same order, Pade approximations are usually superior to Taylor expansions especially when 
functions contain poles. 
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§5. Changes in the image properties 

For a point mass lens model, we have the following "universal" relations for the 
image properties. First, the sum of the two image positions is 

Namely, the sum of the image positions is always equal to the original source position. 
Second, the difference of the amplification factors is always unity: 

Af s = A+-A^ = l. (5-2) 

This also means, for a point mass lens, the image with positive parity is always 
brighter than that with negative parity. Just for reference, we also show that the 
image separation and the total amplification are expressed in the following way: 

Ax = x£-x = y/ ^ +4 /?x, Ayo = Vo ~ Vo = ( 5 ' 3 ) 

B 2 + 2 

4 ot = A+ + An = V . (5-4) 

In this section, we investigate how the quadrupole moment changes the above 
image properties. From Eqs. (|3-9[) and (|3-10p . we obtain the sum as 

x+ + x-=p x -e^ + j^\p x , V + + y- =Py + e(l + ^pjf3 y , (5-5) 

The change due to e is getting larger as oc 4e/3" 2 for \/3\ < 1. Therefore, even if 
|e| <C 1, the "universal" relation Eq. (|5-ip may be significantly broken if <C 1. 



Table III. Comparison of our approximate amplification factors with numerical ones 



e 


& 


Py 


parity 


(l) -<4num 


(2) A 


1(2) -(1)1/(1) 


(3) A P 


1(3) -(1)1/(1) 


0.01 


0.0 


0.2 


+ 


2.83741 


2.82454 


0.45% 


2.83847 


0.04% 










2.37779 


2.32454 


2.24% 


2.37166 


0.26% 


0.01 


0.0 


0.5 


+ 


1.56603 


1.56568 


0.02% 


1.56609 


0.00% 










0.65303 


0.64568 


1.13% 


0.65116 


0.29% 


0.02 


0.0 


0.2 


+ 


2.66004 


2.61173 


1.82% 


2.66404 


0.15% 










2.87313 


2.61173 


9.10% 


2.83724 


1.25% 


0.02 


0.0 


0.5 


+ 


1.54134 


1.53994 


0.09% 


1.54156 


0.01% 










0.73395 


0.69994 


4.63% 


0.72434 


1.31% 


0.01 


0.2 


0.0 


+ 


3.26482 


3.25015 


0.45% 


3.26619 


0.04% 










1.78910 


1.75015 


2.18% 


1.78564 


0.19% 


0.01 


0.5 


0.0 


+ 


1.61749 


1.61714 


0.02% 


1.61758 


0.00% 










0.54292 


0.53714 


1.06% 


0.54170 


0.22% 


0.02 


0.2 


0.0 


+ 


3.52534 


3.46296 


1.77% 


3.53232 


0.20% 










1.59981 


1.46296 


8.55% 


1.58928 


0.66% 


0.02 


0.5 


0.0 


+ 


1.64434 


1.64288 


0.09% 


1.64460 


0.02% 










0.50371 


0.48288 


4.14% 


0.49971 


0.80% 
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On the contrary, the image separation is 

Ax = x + -x- = Ax < 1 - ( 7o'7 o2 ' V A [ ~ 1 I e J> , (5-6) 




4 »=» + -»" =/1 *\ 1 -^^+5 +1 j7- ( "> 

It is evident that, even if <C 1, the change due to e in the image separation is 
~ e x O(l). Therefore, we may conclude that the change due to e in the image 
separation is small for |e| <C 1. 

Similar situation occurs in the analysis of the image amplification. From Eqs. (|4-5jl 
and (|4-6|) . the amplification difference for |/3| <C 1 is 

A dlS = A+-A p ~ 2 (5-8) 

1 + le' - (§ e' * 

where 

'-((*)'-(*)> 

The change due to e is getting larger as oc 2e/3 -2 for |/3| < 1. Therefore, even if 
|e| 1, the "universal" relation Eq. ()5 2[) may be significantly broken if \j3\ ^ 1. 
As is already shown numerically in Table IIH[ there is even a case of ^4 dlff < for 
e = 0.02, 13 X = 0.0, p y = 0.2. 
The total amplification is 



/(tot 

A tot = Aj, + A p ^ °— -J. (5-10) 

2 . ' " 



1 + ¥ ~ {ft 

Compared to ^4 dlff , the change due to e is relatively mild for |/3| <C 1. 

We give some comments on the range of the parameters e and (3. Since we 
have started from the linear perturbation with respect to e, the range of validity is 
< e <C 1. If we assume additionally that the maximal error in the amplification 
factors should be less than, say, 1%, the range of validity may be < e < 0.01, 
which we can guess from Table IIIH The ranges of (3 X and (3 y are more important. 
From Eqs. (|5-8|) and (|5-10|) . we can obtain the constraint |e'| < f3 2 and |e'| <C /3. If 
not, then the term of 0(e 2 ) in the denominators of Eqs. (|5-8p and (|5T0p can make a 
dominant contribution over the linear term, which is outside of the validity of linear 
perturbation theory. 

§6. Summary 



In this paper, we have investigated the gravitational lens effect caused by a non- 
spherical, compact object. The non-spherical property of the gravitational potential 
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is expressed by the quadrupole moment. Eqs. (|2-6p and (j2-7|) are the lens equations 
for a compact lens object with a quadrupole moment. In these lens equations, 
the effect of the non-spherical contribution is expressed by e, the eigenvalue of the 
normalized quadrupole moment tensor Qij. 

Under the assumption that the quadrupole moment e is small, we have pertur- 
batively solved the lens equation and obtained the solutions for the image positions. 
Eqs. (|3-9p and (|3-10p represent the perturbative solutions for the positions of the 
"major" images, which are reduced to the solutions for a point mass lens model in 
the limit e — > 0. We have also found two new perturbative solutions for the "minor" 
images, which are given in Eq. (|3- 18j) . The minor images are very dim, always appear 
very close to the y-axis, near (0, zt-^/e), and vanish in the limit e — > 0. The accuracy 
of our approximate solutions is numerically checked and summarized in Tables U and 
IITI for some sets of parameters. It is shown that the maximal error in the image 
position is about 10 -4 , which is the same order of magnitude as 0(e 2 ). 

Then, we have calculated the amplification factors for the major images from 
the Jacobian and obtained the Pade approximants of them in Eqs. (j4-5|) and ()4-6p . 
The accuracy of our approximate formula for the amplification factor is numerically 
estimated and summarized in Table ITO1 for some sets of parameters. It is shown that 
the typical relative error is 1% or less. 

For a point mass lens model, we know the simple "universal" relations for the 
image properties, Eqs. (|5-ip and (|5 2|) . We have investigated how the quadrupole 
contribution e changes such "universal" relations. We have found that the change 
due to a non-zero e in the image separation is ~ e x 0(1). However, the difference 
of the amplifications A dlS = A + — A~, which is always unity in the case of a point 
mass lens, may be significantly changed due to the quadrupole contribution, which 
is shown in Eq. ()5-8|) . 

Mao and Schneider^ discussed the "anomalous" flux ratio between the images 
caused by a lens galaxy. The lens system they discussed was QSO B 1422 + 231. 
Whereas the well-known "simple" lens models, such as a singular isothermal ellipsoid, 
could fit the observed image positions very accurately, they all failed to obtain the 
observed flux ratios. The possibility they argued is that the discrepancy between the 
observed and model-predicted flux ratios is due to substructure or perturbation in the 
simple lens models, which only slightly changes the image positions but significantly 
affects the flux ratios. Since the lens model in this paper is different from that of 
Mao and SchneiderpSJ we have no intention of directly applying our results to their 
lens systems. Still, it should be noted that the results in our paper show some formal 
similarity. It indicates that even a tiny non-spherical distortion of the lens potential 
may also cause significant amount of flux anomalies in the lensed images, whereas it 
only slightly changes the image positions. 
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